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Active galactic nuclei and gamma ray bursts at cosmolog- 
ical distances are sources of high-energy electron and muon 
neutrinos and provide a unique test bench for neutrino insta- 
bility. The typical lifetime-to-mass ratio one can reach there 
is r/m ~ 500 Mpc/ cE v ~ 500 s/eV. We study the rapid decay 
channel Vi — ► Vj + 4>, where <j> is a massless or very light scalar 
(possibly a Goldstone boson), and point out that one can test 
the coupling strength of gtjViVj down to gtj < 10 -8 eV/m by 
measuring the relative fluxes of i/ e , and v T . This is orders 
of magnitude more stringent bound than what one can obtain 
in other phenomena, e.g. in neutrinoless double beta decay 
with scalar emission. 
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By measuring the spectra of the neutrinos pro- 
duced in Earth's atmosphere by cosmic rays the Super- 
Kamiokande experiment has recently got evidence that 
the muon neutrino has a mass and that it mixes with 
another neutrino flavour [jjj. The observed deficit and 
angular dependence of the atmospheric flux can be 
explained in terms of — v T or — v s oscillations 
if the squared mass difference of the oscillating states 
is Am 2 = 3 • 1(T 4 - 8 • 10~ 3 eV 2 and their mixing is 
sin 2 26 > 0.8. Here v s is a sterile neutrino which lacks 
the standard electroweak interactions. 

One may draw from the Super-Kamiokande result the 
conclusion that in general neutrinos have a mass and lep- 
ton flavour numbers are not separately conserved. This 
would consequently imply that not all neutrinos are sta- 
ble particles but may decay into lighter neutrinos and 
other light particles. If the oscillating partner of is 
v T , the corresponding mass eigenstates are (ignoring all 
other possible mixings) v\ — cosO — sin# v r and 
V2 = s'm9 + cos 6 v T , where 32° < 9 < 58°, and as- 
suming 77i2 > mi ' then the radiative decay z/2 — ► Vi +7 is 
possible. Given the small mass difference of the neutrinos 
as indicated by the Super-Kamiokande result, this decay 
is in general very slow as it proceeds through loops and is 
suppressed due to GIM mechanism and helicity match- 
ing. Assuming the Standard Model (SM) interactions 
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one goes beyond the SM, other decay modes with shorter 
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lifetimes may appear Q . The shortest lifetimes typically 
arise in majoron models B , for decays with a light neu- 
trino and a massless scalar as decay products. 

The question now arises how one could detect the pos- 
sible decays of neutrinos and thereby obtain independent 
new evidence of their masses and mixings and informa- 
tion on their interactions. In this letter we shall address 
this question by considering the decays 
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where vi and vi are light and heavy mass eigenstate neu- 
trinos (or antineutrinos) and <fi is a massless or very light 
scalar. In the simplest majoron model <f> is a massless 
Goldstone boson, a majoron, associated with the spon- 
taneous breaking of global lepton number symmetry. In 
the following we will not restrict ourselves to this specific 
model but will consider more general forms of interac- 
tions where the coupling strength may be independent 
of the origin of the neutrino mass. Our main result is 
that the decay (0) could be detected by measuring the 
flavour content of neutrino flux from active galactic nu- 
clei (AGN) or gamma-ray bursts (GRB's). We will find 
that this will provide several orders of magnitude more 
sensitive probe of the coupling strength than other ex- 
periments. 

Let the scalar-neutrino interaction to be of the form 



>C in t = g l jV c i Vj(l> + h.c. 
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As the scalar <j> should be very light, it has to be dom- 
inantly a weak isospin singlet in order to satisfy the 
constraint from the LEP measurement of the invisible 
Z width ||. If the lepton number L is conserved, the 
scalar has L(<p) = —2. Examples of this kind of mod- 
els are the so-called charged majoron models H and the 
model where Dirac neutrinos couple to a scalar through 
their singlet right-handed components via the coupling 
v\jVjji<f> 0| . In these models the coupling strength gij 
is not related to neutrino masses and there are no con- 
straints from the neutrinoless double beta decay as the 
lepton number is conserved. If the lepton number is bro- 
ken, so that L{4>) — 0, neutrinos are in general Majo- 
rana particles allowing for neutrinoless double beta de- 
cay. This is the case in the majoron model. 

In the simplest version of the model the coupling g 
is proportional to neutrino mass, and the ensuing decay 
time is much too long to be relevant for our discussion. 
However, other models can be easily constructed where 
the interaction matrix is not proportional to the mass 
matrix and where shorter lifetimes are possible S. 
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Taking 500 Mpc as a typical distance (D) of AGN or 
GRB's and 500 TeV as a typical neutrino energy (E), the 
typical lifetime-to-mass ratio one can reach when consid- 
ering decays of AGN or GRB neutrinos is 
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which is valid for all experimentally allowed neutrino 
mass values. 

The width of the decay Vj — > v f + in the rest frame 
of v,j is given by 
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where <5to 2 ,- = mj — to 2 . The lifetime in the rest system 
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of the observer, where Vj has the energy Ej, is then 
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For a relativistic decaying neutrino the decay distance is 
then given by 
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In order to see the effect of the decay, one should ob- 
serve a sufficiently abundant neutrino beam from a flight 
distance of D > L. As the numbers in (^J) indicate, as- 
trophysical or cosmological distances are preferred, and 
AGN and GRB's stand out as possible beam sources. 
Note, however, that neutrino decays may be relevant also 
with respect to the atmospheric neutrino problem where 
smaller distances and energies are involved M . The high- 
est energy cosmic rays are known to originate in AGN 
and GRB's, and they are assumed to produce also neu- 
trinos from the decay of charged pions through the chain 



r + 



,+ 



cnarged 



7r — > is^/i — > v^e ' Vf-v^ |[1U|| . The primary process 
where the pions are produced is the interactions of ac- 
celerated protons with the photons via p7 — * mt + in the 
source near the A resonance. The energies of neutri- 
nos originating in AGN jets are typically 10 5 TeV, but 
less energetic neutrinos are abundantly produced in ac- 
cretion disks |Oj. The typical energies of GRB neutri- 
nos are around 100 TeV. The most powerful AGN and 
GRB's are located at distances of 100-1000 Mpc or more. 
AGN and GRB neutrinos whose energy exceeds 10 TeV 
can be distinguished from the atmospheric neutrino back- 
ground jL2}. 

If we take E = 100 TeV, D = 500 Mpc as face val- 
ues, we will find from (^) that the values of the coupling 
constant down to 
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could in principle be probed by observing the neutrino 
flux from AGN and GRB's and measuring the relative 
fractions of its flavour components. 

The neutrino beam created at AGN and GRB's con- 
sists of v e and with no substantial v T component. If 
we assume that the oscillating partner of atmospheric 
is the tau neutrino, oscillations will create a v r compo- 
nent into the AGN and GRB neutrino beam of about the 
same size as the component as the mixing is almost 
maximal. Actually, with such a baseline the presence of 
a v T component would indicate neutrino oscillations with 
Am 2 as low as 10~ 17 eV 2 HQ. 

In the case of almost maximal mixing the decay will 
not dramatically change the content of the beam from the 
fifty-fifty division as both decaying and produced mass 
eigenstate neutrinos contain and v T components in 
almost equal fractions. At most, the abundance ratio of 
Vrlv^ in the beam is changed due to the decay by the 
factor cot 9 ~ 1.6, corresponding to the lowest mixing 
angle allowed by the Super-Kamiokande measurement. 

In contrast, a very clear signal would be available in the 
case where the decay happens into a third neutrino mass 
eigenstate consisting mostly of the electron neutrino v e 
but having no substantial and v r components. Actu- 
ally, in the case where the — v r mixing is responsible 
for the atmospheric neutrino signal, the electron neutrino 
should mix with a sterile neutrino in order to explain the 
solar neutrino deficit problem, i.e. there should exist an- 
other two mass eigenstates ^3 = cos <j) v e — sin cf> v s and 
z/4 = sin0 v e + cos <p Vsi where v s is a sterile neutrino. 
Three solutions are still possible; two of them are based 
on the MSW mechanism Jl5[ | , one with small mixing an- 
gle (sin 2 29 = 5.5 x 10~ 3 ), one with large mixing angle 
(sin 2 29 ~ .8) , and the third one is based on vacuum 
oscillations [[l6| with large mixing angle (sin 2 29 ~ 0.75). 
For the MSW solutions the squared mass difference Am 2 4 
is 5 x 10~ 6 eV 2 and for the vacuum oscillation solution 
10~ 10 eV 2 . The most natural mass hierarchy would be 
that mi ~ to 2 3> to 3 ~ TO4. Because of the strong mixing 
of and v T , both heavy components of the original AGN 
neutrino beam, v\ and V2, would in this case be expected 
to decay into the lighter mass states V3 and 1^4, resulting 
in the disappearance of the expected AGN muon neutri- 
nos, as well as the non-appearance of tau neutrinos in 
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If the atmospheric neutrino behaviour is due to mixing 
between and a sterile neutrino u s , then the solar neu- 
trino deficit problem could be solved through a v e — v T 
mixing. This case differs from the previous one in that 
now the disappearance of muon neutrinos could be ac- 
companied with the appearance of tau neutrinos as a 
result of decays. If the v e — v T mixing is large, as is the 
case when the solar neutrino deficit results from vacuum 
oscillations or a large- mixing MSW effect, then a large v T 
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component would be generated also through v e — v r os- 
cillations, but in the case of the small- mixing MSW effect 
z/j-'s would appear almost exclusively due to decays. 

To experimentally observe the possible decays of AGN 
and GRB neutrinos it is hence essential to measure the 
relative fluxes of v e , and v T . While the muon neutrino 
will be relatively easy to detect by tracking the muon 
produced in charged-current interaction in the vicinity 
of a Cherenkov detector in deep water or ice, the de- 
tection of electron and in particularly tau neutrinos is 
less straightforward. The electron neutrinos create an 
electromagnetic cascade which can be detected by op- 
tical or radio technique W\ - For the detection of the 
tau neutrino two methods has been recently proposed: 
by double-bang events from the production and decay 
of the r lepton fl3[| , and by the absence of absorption 
by the Earth (TJJ] . In general, the effective elimination of 
the background requires that one should look at up-going 
neutrinos. 

Depending on the scenario realized in Nature, one 
can get information on various coupling constants gij 
(i,j = 1,2,3,0) or their combinations by measuring the 
AGN and GRB neutrino flux in the Earth. As men- 
tioned already, the sensitivity is in principle on the level 
of < 10~ 8 . Limits obtained from other astropysical 
objects, such as supernovas, and from early Universe 
are in general less stringent (see Ref. [|l8| for a review). 
The strongest laboratory bounds on the scalar couplings 
9ij v i v j are due to the scalar-emitting neutrinoless dou- 
ble beta decay They apply to the couplings where 
both the mass eigenstates Vi and Vj contain v e flavour 
component, and the magnitude of the upper bound for 
gij depends on whether the scalar is a Goldstone boson 
or not as this affects the form of the electron spectrum. 
In the case 4> is not a Goldstone boson the constraint 
is g^ < 10~ 4 and in the case it is g^ < 0.1. These lim- 
its do not compete with the limits one can achieve from 
AGN and GRB's. Let us note that exotic neutrino in- 
teractions of the AGN and GRB neutrinos can be also 
probed via their interactions with the cosmic neutrino 
background [^o| . 

To summarize, we have pointed out that if the neu- 
trino beam from AGN and GRB's contains neither v^s 
nor ^ r 's the only explanation in the framework we have 
considered would be the decay of the containing mass 
eigenstates into a lighter neutrino and a scalar. If one 
observes a v r component in the neutrino flux but no 
component this will also indicate neutrino decay. If there 
is also a component in the beam, then the v T compo- 
nent could also result from oscillations. Because of large 
distance to the source this test of neutrino stability is or- 
ders of magnitude more sensitive than other astrophysical 
and cosmological, as well as laboratory, tests. 
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